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ABSTRACT 

Using data obtained for twelve galaxies as part of the Spitzer Infrared Nearby Galaxies Survey (SINGS) 
and the Westerbork Synthesis Radio Telescope (WSRT)-SINGS radio continuum survey, we study how star 
formation activity affects the far-infrared (FIR)-radio correlation within galaxies by testing a phenomenolog- 
ical model which describes the radio image as a smeared version of the FIR image. The physical basis of 
this description is that cosmic -ray (CR) electrons will diffuse measurably farther than the mean free path of 
dust-heating photons before decaying by synchrotron radiation. This description works well in general. Galax- 
ies with higher infrared surface brightnesses have best-fit smoothing scale-lengths of a few hundred parsecs, 
substantially shorter than those for lower surface brightness galaxies. We interpret this result to suggest that 
galaxies with higher disk averaged star formation rates have had a recent episode of enhanced star formation 
and are characterized by a higher fraction of young CR electrons that have traveled only a few hundred parsecs 
from their acceleration sites in supernova remnants compared to galaxies with lower star formation activity. 

Subject headings: cosmic rays — infrared: galaxies — radio continuum: galaxies 
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1. EVITRODUCTION 

Radio continuum emission from normal galaxies arises 
from a combination of non-thermal (synchrotron) and ther- 
mal (free-free) processes; the former critically depends on a 
galaxy's magnetic field and cosmic ray (CR) electron distribu- 
tions. To date, the evolution of CRs has been studied directly 
only in the Milky Way. In other galaxies, information has 
come from sampling the CR electron component via multi- 
frequency radio observations (e.g. Duric 1991 ; Dahlem et al. 
fT^gSHLisenfeld et al. 1996a; Irwin et al. 1199 9: Beck 2005). 
Radio continuum data alone, however, only describe the 
present distribution of CR electrons, without providing any 
information regarding the initial source distribution or propa- 
gation history. 

The close spatial correlation between the thermal dust 
and non-thermal synchrotron emission within galaxies ar- 
gues for a shared origin presumably massive star formation 
iHarwit & Pacini I Il975t iHelou. Soifer. & Rowan-Robinson I 
Il985h . Young massive stars, still in or near their dusty na- 
tal molecular cloud, emit photons that are re-radiated by dust 
grains at far-infrared (FIR) wavelengths. These stars are also 
the progenitors of supernovae (SNe) whose remnants (SNRs) 
are thought to accelerate CR electrons responsible for syn- 
chrotron radiation. Since the mean free path of dust-heating 
photons (~100 pc) is far shorter than the diffusion length of 
CR electrons (^1-2 kpc), Bica v & Jie lou ( 1990) conjectured 
that the radio image of a galaxy should resemble a smoothed 
version of its infrared image. This prescription has been 
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shown to hold for galaxies observe d at the " super resolution" 
(<1') of IRAS HIRES data (Marsh & Helou 1998). Using 
Spitzer imaging, Murphv et al. (2006, hereafter M06) tested 
this model on sub-kiloparsec scales within four nearby galax- 
ies and found it applicable to first order, improving the corre- 
lation between the infrared and radio maps while leaving pat- 
terns in the residuals th at mirror the st ructure in the galaxy's 
image. More recently, iHughes et al. I li2006) has found that 
within the Large Magellanic Cloud (LMC), synchrotron ha- 
los around individual star-forming regions are more extended 
than FIR-emitting regions, corroborating this phenomenology 
on scales >50 pc. 

In this Letter we explore how star formation activity affects 
the relative spatial distributions of FIR and radio continuum 
emission for a sample of 12 spiral galaxies observed as part of 
the Spitzer Infrar ed Nearby Galaxies Su rvey (SINGS) legacy 
science project (Kennicutt et al. '2003*) and the Westerbork 
Synthesis Radio Telescope (WSRT)-SINGS radio continuum 
survey. This analysis builds on that of M06, though the larger 
sample has increased the range in star formation rates (SFRs) 
investigated by a factor of Such a study is only now pos- 
sible due to the high spatial resolution and sensitivity of the 
Spitzer Space Telescope at FIR wavelengths. 

2. DATA AND ANALYSIS 
2. 1. Observations and Data Reduction 

Observations at 24, 70, and 160 /jm were obtained us- 
ing the Multiband Imaging Photometer for Spitzer (MIPS; 
iRieke. et al. 2004) as part of the SINGS legacy science pro- 
gram. The MIPS data were processed using the MIPS Data 
Analysis Tool (DAT; Gordon et al. 2005) and included in the 
SINGS data release 3 (DR3). Details on the data reduction can 
be found in M06. The full width at half maximum (FWHM) 
of the MIPS 24, 70, and 160 ^m point spread functions (PSFs) 
are 5"7, 17", and 38", respectively, and the final calibration 
uncertainty is ^10 % at 24 /jm and ~20 % at 70 and 160 /jm. 

Radio continuum imaging at 22 cm was performed us- 
ing the WSRT. Each galaxy was observed in a configura- 
tion with particularly good sampling at short baselines to cur- 
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tail flux loss at low spatial frequencies. The WSRT data 
were CLEANed and self-calibrated using an imaging pipeline 
based on the MIRIAD package. The final CLEAN maps were 
restored with 18" FWHM circular Gaussian beams. Total 
intensity calibration of the data was performed in the AIPS 
package and the flux density calibration of the radio maps is 
better than 5 %. A c omplete description of the image process- 
ing can be found in lBraun et al. I(l2006h . 

To properly match the resolution of the MIPS and radio im- 
ages, the Fourier Transform (FT) of the 70 /jm images was di- 
vided by the FT of a model of the MIPS 70 /jm PSF and then 
multiplied by the FT of a Gaussian PSF matching that used as 
the CLEAN restoring beam. The final product was then trans- 
formed back into the image plane, and checked to ensure flux 
was conserved. By matching the resolution of the MIPS and 
radio data with Gaussian beams, rather than the MIPS 70 /jm 
PSF (which has significant power in its wings), the dynamic 
range of the 70 /jm images has increased by roughly a factor 
~2. Due to the improved resolution there is also a slight in- 
crease in the RMS noise of the maps by an average of ^33 %. 

2.2. Image-Smearing Model 

We smooth each galaxy's 70 fjm map with a parameter- 
ized kernel and calculate the normalized squared residuals be- 
tween the smoothed 70 /jm and the 22 cm image as defined by. 



(1) 



where Q = L^y/L^; is used as a normaUzation factor, I{1) 
represents the infrared image after smearing with kernel of 
scale-length /, R is the observed radio image, and the subscript 
j indexes each pixel. The smoothing kernel takes the form, 
K(r) = e^''!''", where ro is the projected, e-folding length of 
the kernel (1). This kernel was chosen because it works as 
well as or better than the other kernel types studied by M06. 

We calculate the residuals after first removing pixels de- 
tected below the 3 O level in either the radio or maximally- 
smeared infrared maps and editing out contaminating back- 
ground radio sources. Additional editing was performed in 
the cases of NGC 3031 and NGC 5194 as described in M06. 

The residuals are plotted against the smearing kernel scale- 
length in Figure ^ To quantify the improvement in the spa- 
tial correlation between the 70 /jm and 22 cm maps due to 
smoothing the 70 /jm image, we measure <I> defined as the log- 
arithmic decrement from (|)(2,0) to (j)(2iO at its minimum; 
this minimum point defines the best-fit scale length. Since 
<I> is an estimate of the global improvement, we also created 
residual maps, defined by log (2^^/)— log using the best- 
fit smearing kernels to determine how well the smearing pre- 
scription works for different regions within each galaxy. The 
70 jum and residual maps are plotted in Figure|2l 

2.3. SFR Surface Densities 

To see how the FIR-radio correlation within galaxies is af- 
fected by ongoing star formation activity, we measure their to- 
tal infrared (TIR; 3-1100 /jm) surface brightness and express 
it as SFR surface density, Esfr- We first create /tir maps via a 
weighted combination of the 3 MIPS maps according to Equa- 
tion 4 of Dale & Helou ( 2002). The SFR surface densities, in 
units Mp yr^' kpc^^, were then calculated by estimating a 
SFR (i.e. Equation 5 of lBellll2 003) using the TIR luminosity 



within an ellipse fit to the 1.4 x 10^ W mT sr^ isophotal 
radius, and then dividing by the corresponding deprojected 
area (comparable to the galaxy sizes at 70 /jm). 
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Fig. 1. — Residuals betw een the observed radio maps and smeared 70 ^im 
images (as defined in i|2.2l as a function of smearing scale-length. Galaxies 
havin g low disk-averaged star formation rates, defined by log ZsFR < —2.3 
(see j|2.3l are plotted in panel (a) while galaxies with high star formation 
activity, and larger values of ZsFR, are plotted in panel (b). 



3. RESULTS 

In Figure we plot the residuals against the smearing 
scale-length for the sample galaxies, separated into low and 
high Esfr bins. There is a wide range of behavior in the 
residual curves among the galaxies. Three galaxies do not 
display any definitive minimum in their curves. We also find 
that galaxies with higher Esfr have better defined single min- 
ima compared to galaxies with lower Esfr even though the 
distribution of <I> is similar between the two groups. 

The smearing technique improves the correlation between 
the radio and 70 /jm maps by an average of ^0.5 in log (|). 
This value is 0.3 dex larger than what was found by M06 and 
is due to improvements in the resolution matching and the 
larger sample. 

We display the 70 /jm and residual maps of each galaxy 
ordered by Esfr in Figure |2l While there appears to be a 
general correlation between the structure found in the 70 /jm 
and residual maps, there is a range in the residual map be- 
havior exemplified by the cases of NGC 628, NGC 5055, and 
NGC 6946. For NGC 5055, the residuals are small and the 
image is nearly uniform. In NGC 6946, on the other hand, we 
see distinct structures in the residual maps associated with the 
spiral-arm morphology of the galaxy, as well as a much larger 
range in the residuals compared to NGC 5055. The arms of 
NGC 6946 display infrared excesses peaked on star-forming 
regions while the galaxy's inter-arm regions exhibit excess ra- 
dio emission. NGC 628 is distinctly opposite to NGC 6946 
in its residual map displaying radio excesses for star-forming 
regions. This implies that the large smoothing scale-lengths. 
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which provide the best-fits for low Esfr galaxies, oversmooth 
the active star forming components. 

In Figure|3lwe plot the best-fit scale-length for each galaxy 
versus their SFR surface density, Esfr- Galaxies are bet- 
ter modeled by larger smearing scale-lengths if their disk- 
averaged SFR is lower To ensure that this trend is not an arti- 
fact, we have compared the best-fit scale-lengths with galaxy 
distances, linear sizes, and signal-to-noise (S/N) ratios and 
find a correlation with none of them. The implications of this 
trend will be discussed in ^ 

3.1. Effects of Free-Free Emission 

Since the image-smearing technique is designed for the 
physics of non-thermal radio emission, the presence of free- 
free emission can complicate the interpretation of the best-fit 
smearing scale-lengths. To gauge the effect of free-free emis- 
sion on the trend in Figure|3] we repeat the smearing analysis 
after subtracting an estimate of the free-free emission distri- 
bution using a scaled version of the 24 /jm maps. 

The scaling factor was determined by first relating the 
24 /jm emission to Paa line emission using the correlation 
found between 24 /jm and Paa luminosities within NGC 5 194 
(M51a) (Calzetti et al. 2005). Although this correlation is not 
universal ( Perez-Gonzalez et al. „2 QQ( ^,). it should be sufficient 
for the purpose of roughly estimating how the thermal emis- 
sion affects the best-fits scale-lengths. The corresponding Paa 
emissio n^stimate can then be related to an ionizing photon 
rate (Osterbrock 1989), and the expected free-free emission 
at 1.4 GHz (.Rubin -1968) . The final relation is given by, 

\]y ) VlO'^Ky VGHz/ V Jy / 

(2) 

where we have assumed an average H II region temperature 
ofr = 10'^K. 

For seven of the sample galaxies, 'Ni klas et al. I (Il997h re- 
port a mean thermal fraction of ^0.08 at 1.0 GHz from ra- 
dio spectral index fitting. Assuming a mean spectral index 
of —0.8, we obtain a very similar thermal fraction using our 
data. Since the distribution of free-free emission is clumpy, 
mostly originating in H II regions, we expect that the smearing 
scale-lengths should increase after it is removed. The best-fit 
scale-lengths increase by an average of ^25 %. The slope of 
the fit using these scale-lengths versus Esfr is shallower by 
<10 %, indicating that the observed trend in Figure |3l cannot 
be suppressed by removing the thermal radio emission. 

4. DISCUSSION 

We have analyzed the effect of star formation intensity on 
the FIR-radio correlation within galaxies using a phenomeno- 
logical image-smearing model. We find that the best-fit scale- 
length decreases as the star formation activity within each 
galaxy increases, implying a systema tic decrease in the mea n 
distance diffused by CR electrons JHelou & Bicav I ri993l) . 
The measured diffusion scale-lengths can provide indepen- 
dent estimates for the ages of related star formation episodes 
as well as help constrain interstellar medium (ISM) param- 
eters associated with the propagation characteristics of CR 
electrons. 

The physical interpretation of our results indicates that CR 
electrons are, on average, closer to their place of origin in 
galaxies having higher star formation activity. Since the dif- 
fusion scale-length of CR electrons depends only on their age 
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Fig. 3. — The SFR surface densities, (ZsFR. see i|2.3l . for each galaxy 
are plotted against the best-fit smoothing kernel scale-length. The general 
improve ment in the spatial correlation between the 70 fjm and 22 cm maps 
(<3>, see jl2.2t is indicated by the plotting symbol. Also plotted is the fit to 
the data (dotted line) along with the expected diffusion scale-lengths due to 
inverse Compton (IC) losses in a fixed magnetic field (dashed line) and syn- 
chrotron + IC losses with an energy-dependent diffusion coefficient De for 
the steepest possible index (solid line) (see 0. The increase to the best-fit 
scale-length using 22 cm maps correc ted f or free-free emission is indicated 
by a vertical line for each galaxy (see i|3.1l . 

and ability to diffuse through a galaxy's ISM, three explana- 
tions of the results are possible. The CR electrons may: (1) 
have relatively short lifetimes due to a high energy loss rate; 
(2) diffuse at a slower rate due to the ISM having a high den- 
sity and magnetic field strength, resulting in a shorter mean 
free path; or (3) have been accelerated recently and be rela- 
tively young. The first two of these explanations are applica- 
ble in the case of steady-state star formation, while the third 
requires a recent episode of enhanced star formation. We will 
try to distinguish among these three explanations. 

As CR electrons propagate through the ISM, they lose en- 
ergy due to synchrotron, inverse-Compton (IC) scattering, 
bremsstrahlung, ionization, and adiabatic expansion losses. In 
normal star-forming galaxies only the first two mechanisms 
are significant for CR electrons associated with synchrotron 
emission at ^^1 GHz (Condon 1992). Assuming an isotropic 
distribution of electron velocities, ' the cooling timescale due 
to synchrotron and IC losses for CR electrons emitting at a 
frequency v, in a magnetic field of strength B is, 

f ^) ^2.8 X 10^ f ^) f A) ( ^B + ^-a. \ 
V yr y VGHz/ \ijG J V 10"'^ erg cm-V 

(3) 

where Ub and Urad are the magnetic and radiation field energy 
densities, respectively. 

In simple diffusion models, the propagation of CR electrons 
is usually characterized by an empirical energy-dependent 
diffusion coefficient, De, where De = 10^^ cm^ s^^ for 
E < 1 GeV and I02'*(£/GeV)i/2 cm^ s"' for E > 1 GeV 
dGinzburg et al. I '1980). Though derived for CR nuclei, the 
above values appear consistent with empirically measured 
diffusion coefficients for CR electrons (e.g. Dahlem et al. I 
»1995i) . According to synchrotron theory, CR electrons will 
emit at V o<: BE^ thereby allowing De to be expressed as a 
function of B for a fixed v. Now, assuming that CR electrons 
diffuse a distance Zcooi — V^e'^cooI, we can estimate the rel- 

' For isotropically distributed electron velocities we assume < sin"0 >= 
|, where is the CR electron pitch angle. 
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ative importance of IC losses, ISM density (n), and B on the 
distance CR electrons travel. 

Fixing B and scaling f/rad °= J^sfr, we plot the expected re- 
lation between /cool and t/rad in Figure |3] (dashed line). This 
relation has a slope ~9 times shallower than a line fit to the 
data (dotted line). It is more likely that n, B, and will all 
scale together, with a De dependence on n. Accordingly, we 
introduce the scaling relations, 

Boc„P, Z)£oc„5g-i/4^ (4) 

where t < 3 < 1 and — 1 < 5 < — | are realistic index ranges 
(flelou & Bicay^ 1993). 

Setting Ub = U^d, which has been shown to be valid for a 
large sample of spirals cLisenfeld et al 1996b.) . leads to B oc 

^SFR' '^cooi °= B^^l^, and the dependence 

/coo,c<Er/r-'/'''- (5) 

For the unlikely case that each galaxy takes the extreme values 
of P = ^ and 5 = — 1, we obtain the steepest possible index of 
— 19/16 (solid line in Figure |3}, which is still a factor of ^2 
times shallower than the line fit to the data. Consequently, dif- 
ferences in the diffusion properties and/or cooling timescales 
of each galaxy's CR electron population alone (i.e. the case 
of steady-state star formation), cannot explain the observed 
trend in Figure|3] 

On the other hand, by setting the best-fit scale-lengths equal 
to /cool and B = 5 /jG for each galaxy, we estimate mean ages 
ranging from ^1.1x10^ to ^8.3x10^ yr for these 1.4 GHz 
emitting populations of CR electrons. Even if we assume 

1/2 

B °^ J^sFR, this range of ages shrinks only slightly. While these 
ages are plausible, spanning nearly 3 orders of magnitude 
over a factor of 4 change in Esfr seems unphysically large. 
This suggests that as Esfr decreases the best-fit scale-length 
shifts from being dominated by young CR electrons associ- 
ated with regions of star formation to being dominated by old 
CR electrons in the galaxy's underlying disk. Consequently, 
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we conclude that galaxies with higher Esfr are likely to have 
experienced a recent episode of enhanced star formation com- 
pared to galaxies with lower disk-averaged SFRs, leading to 
a larger fraction of relatively young CR electrons within their 
disks. We should also note that galaxies with higher star for- 
mation activity may also suffer from increased CR electron 
escape due t o instability driven breaks in magnetic field lines 
JParkerlll96 6). though characterizing the importance of es- 
cape is beyond the scope of this Letter. 

The galaxies analyzed here exhibit a wide range of ages and 
spreading scale-lengths of their CR electrons, going from the 
youngest populations, recently injected within star-forming 
regions (/ ^^500 pc), to the oldest CR electrons, making up 
a galaxy's underlying synchrotron disk (/ ^3 kpc). This ex- 
ponential synchrotron disk results from the superposition of 
older CR electrons originating from all star-forming regions. 
In spite of an occasionally dominant mode, a given galaxy 
will contain a range of ages for CR electrons, so it is no sur- 
prise that a single smoothing kernel cannot describe perfectly 
the relation between infrared and radio disks, as evidenced by 
the recognizable structure in the residual maps. NGC 628, 
NGC 2403, NGC 3031, and NGC 3184 exhibit inflection 
points in their residual curves which is also suggestive of mul- 
tiple CR electron populations, possibly from time-separated 
bursts of star formation. A more realistic description to ac- 
count for the range of ages translates into a multiscale analy- 
sis of the galaxy images to separate the signatures of various 
CR electron populations. This will be the subject of a forth- 
coming paper. 
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Fig. 2. — The 70 ;um (left) and residual (right) images are plotted for each galaxy and ai'e ordered by increasing ZsFR from left to right, top to bottom. The 70 /.tm 
images are logaiithmically scaled with a stretch ranging from the RMS background level to the maximum surface biightness in the galaxy disk. Streaking in the 
70 fjm images, clearly visible for NGC 4736, NGC 5055, and NGC 7331, is a latent image effect that does not affect the results (see M06 for details). Residual 
maps are created by subtracting the observed radio maps from the smeared 70 fjm images (as defined in i|2.2l for each galaxy using the best-fit exponential kernel 
projected in the plane of its galactic disk. Regions removed for the residual calculations (e.g. background radio sources) appear as dark, circular holes in a few 
of the residual maps. 



